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ABSTRACT 

We report a study to constrain the fraction of baryonic matter in the cold plus hot dark 
matter (CHDM) universe by numerical simulations which include the hydrodynamics 
of baryonic matter as well as the particle dynamics of dark matter. Spatially flat, 
COSi^normalized CHDM models with the fraction of hot component il^ < 0.2 are 
considered. We show that the models with h/n/flh = 0.5/0.9/0.1 and 0.5/0.9/0.2 give 
a linear power spectrum which agrees well with observations. Here, h is the Hubble 
constant in unit of 100 km/s/Mpc and n is the spectral index of the initial power 
spectrum. Then, for the models with h/n/Clh = 0.5/0.9/0.2 and baryonic fraction f2f, = 
0.05 and 0.1 we calculate the properties of X-ray clusters, such as luminosity function, 
temperature distribution function, luminosity-temperature relation, histogram of gas 
to total mass ratio, and change of average temperature with redshift z. Comparison 
with the observed data of X-ray clusters indicates that the model with = 0.05 is 
preferred. The COBE- normalized CHDM model with rij, > 0.1 may be ruled out by 
the present work, since it produces too many X-ray bright clusters. 
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1 INTRODUCTION 

The large-scale structure of the universe forms when matter 
accretes onto high density perturbations via gravitational 
instability. Among the models that attempt to explain the 
quantitative features of the large-scale structure, the cold 
dark matter (CDM) model is based on the assumption that 
the dominant component of matter is cold. The model also 
assumes that the initial fluctuations are adiabatic and ran- 
dom Gaussian, and have a Zel'dovich spectrum, P(k) = Ak n 
with n ~ 1. The standard CDM (SCDM) model, which is the 
simplest and yet the most widely studied, makes the addi- 
tional assumptions that the total Q — 1 and fit ~ 0.05 from 
the standard big bang nucleosynthesis (Walker et al. 1991). 
Here Q is the density parameter. This model had some suc- 
cess, but it is now well known to have difficulties in explain- 
ing a number of observations (see, e.g., Ostriker & Stein- 
hardt 1995). In particular, this model has excessive power 
on small scales when normalized to the COBE data on large 
scales (see, e.g., Bunn, Scott & White 1995). Several alterna- 



t E-mail: ewchoi@vega. kyungpook. ac.kr; 
ryu@sirius. chungnam. ac 



tive cosmological models, including the CHDM model, have 
been studied claiming better agreement with observations. 

The CHDM model is based on the same assumptions 
as the SCDM model, except that in addition to the CDM 
component, it has also the hot dark matter (HDM) compo- 
nent. HDM is composed of one of the three known species of 
neutrinos {v T , u^, u e ) with a mass of 91.5fih/i 2 eV . Since neu- 
trino's free streaming erases fluctuations in density on scales 
from galaxies to clusters during the radiation-dominated era, 
replacing a portion of CDM with HDM suppresses power on 
small scales. Therefore, the CHDM model has less power on 
small scales than the SCDM model. This basic property of 
the CHDM model was investigated some time ago (Fang, Li 
& Xiang 1984; Valdarnini & Bonometto 1985; Achilli, Oc- 
chionero & Scaramella 1985). 

The fact that the CHDM model is promising in ex- 
plaining the observed data of large-scale structure was first 
established by several analytic calculations using linear or 
nonlinear tools (Schaefer, Shafi & Stecker 1989; Van Dalen 
& Schaefer 1992; Holtzman & Primack 1993; Pogosyan & 
Starobinsky 1995; Dodelson, Gates & Stebbins 1996; Ma 
1996; Borgani et al. 1995, 1997). For example, Holtzman & 
Primack (1993) used the peaks formalism for Gaussian den- 
sity fields and found that the correlation function of clusters 
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in the CHDM model with Qh — 0.3 is consistent with the 
correlation function of the Abell clusters. 

A number of numerical simulations further confirmed 
the promising aspects of the CHDM model (Davis, Summers 
& Schlegel 1992; Klypin et al. 1993; Jing et al. 1994; Klypin 
& Rhee 1994; Nolthenius, Klypin & Primack 1994; Jing & 
Fang 1994; Yepes et al. 1994; Bryan et al. 1994b; Ghigna 
et al. 1994; Bonometto et al. 1995; Walter & Klypin 1996; 
Klypin, Nolthenius & Primack 1997; Dave et al. 1997). Most 
of these simulations were carried out for the model with 
the "standard" density ratio Q c /Qh/Qb = 0.6/0.3/0.1. For 
example, Klypin & Rhee (1994) calculated the correlation 
function for the Abell clusters and the APM clusters and 
confirmed the results given by Holtzman & Primack (1993). 
Jing & Fang (1994) calculated the evolution of the mass 
function, the velocity dispersion function, and the tempera- 
ture function of clusters and argued that their results favor 
the CHDM model over others. Bryan et al. (1994b) com- 
puted the properties of X-ray clusters using a code which can 
handle the hydrodynamics of baryonic matter. They showed 
that the luminosity and temperature distribution functions 
fit well the available observational data, and Qb ?S 0.1 is suf- 
ficient to explain the X-ray properties in the CHDM model 
with Q h = 0.3. Klypin et al. (1997) showed that the CHDM 
model with the "standard" density ratio, as well as that with 
Qh = 0.2, give good fits to a wide variety of "present-epoch" 
data. 

However, the CHDM model with Q. c /Q.h/£lb = 
0.6/0.3/0.1 has a problem in producing galactic halos at 
z <; 3 massive enough to account for neutral gas observed in 
damped Lya systems (Mo & Miralda-Escude 1994; Kauff- 
mann & Chariot 1994; Ma & Bertschinger 1994; Klypin et 
al. 1995; Bi, Ge & Fang 1995). So it was suggested that the 
more promising version of the CHDM model should have 
less mass in the hot component, Q h <J 0.2, in order to have 
more power on galaxy scales and so produce significantly 
more high-redshift objects. Ma (1996), using an approxi- 
mation to the evolution of linear power spectrum, argued 
that the slightly tilted (n ~ 0.9 - 0.95) CHDM model with 
Qh ~ 0.1 — 0.2 gives a power spectrum that agrees best with 
the observed power spectrum. Primack et al. (1995) showed 
that dividing the neutrino mass between two species of neu- 
trinos [v T , i/ M ) lowers the power spectrum on the cluster 
scale, and thus lowers the cluster abundance without the 
necessity of tilt n < 1 of the CHDM model. 

In this paper, we constrain Qb in the CHDM models 
with Qh < 0.2 by investigating the properties of X-ray clus- 
ters that are sensitive to Q b . Recent satellite X-ray observa- 
tions made the properties of X-ray clusters increasingly im- 
portant to cosmology as a probe into the large-scale struc- 
ture of the universe. Being massive and rare, the cluster 
abundance in the local and distant universe carries vital in- 
formation on the initial density fluctuations and the matter 
content of the universe. Also being relatively young dynam- 
ically, the details of their structures provide us with some 
signatures left over from the formation epoch as well as in- 
formation on the background cosmology. The study on the 
properties of X-ray clusters in model universe has been made 
analytically (e.g., Kitayama & Suto 1996) or numerically us- 
ing grid-based codes (Kang et al. 1994; Cen & Ostriker 1994; 
Bryan et al. 1994a; Bryan et al. 1994b) and the SPH code 
(e.g., Evrard, Metzler & Navarro 1996). A general consen- 



Table 1. Model Parameters. 



Parameter 


CHDMa 


CHDMb 


CHDMc 


CHDMd 


h 


0.5 


0.5 


0.5 


0.5 


n 


0.9 


0.9 


0.9 


0.9 


n c 


0.75 


0.7 


0.85 


0.8 


n h 


0.2 


0.2 


0.1 


0.1 


n b 


0.05 


0.1 


0.05 


0.1 


era 


0.67688 


0.67688 


0.73921 


0.73921 



sus seems to be that the SCDM model with h — 0.5 nor- 
malized to the CO BE measurement of the anisotropics in 
the cosmic background radiation (i.e., a$ > 1) has serious 
difficulties in explaining the observed properties of X-ray 
clusters, such as the cluster abundance (Kang et al. 1994; 
Bryan et al. 1994a), the baryon fraction in clusters (Lubin 
et al. 1995), and the contribution of cluster emission to the 
X-ray background (Kang et al. 1994; Bryan et al. 1994a; 
Kitayama & Suto 1996). On the other hand, a flat CDM 
model with cosmological constant (ACDM) and a CHDM 
model, both of which have smaller values of as, seem more 
consistent with observations (Cen & Ostriker 1994; Bryan 
et al. 1994b). Here, ag is the present-epoch, linear rms mass 
fluctuation in the sphere of radius 8/i _1 Mpc. 

We should mention that our work bears many similar- 
ities with that of Bryan et al. (1994b). Both consider X- 
ray clusters in the CHDM universe with similar (but not 
same) codes. However, while Bryan et al. (1994b) considered 
the CHDM model with Q c /Q h /Q b = 0.6/0.3/0.1, which was 
"standard" at that time, we considered the models with a 
smaller Qh {Qh = 0.2, see §3.2), which are more acceptable 
now. 

In §2 the details of the numerical simulations are de- 
scribed. Results are presented in §3. Conclusion follows in 
§4. 



2 NUMERICAL SIMULATIONS 

2.1 Choice of Models 

In order to study the linear power spectrum using analytic 
approximations in the CHDM universe, Ma (1996) consid- 
ered the models with 0.5 < h < 0.8, 0.8 < n < 1, and 0.05 < 
Qh < 0.3. The normalization of the power spectrum, as, was 
determined by the COBE observation (see §2.2 for details). 
As mentioned in Introduction, Ma (1996) found that the 
models with n ~ 0.9 — 0.95 and Q h ~ 0.1 — 0.2 give power 
spectra that agree well with observation. Here, we adopt 
Ma's models with h/n/Q h = 0.5/0.9/0.1 and 0.5/0.9/0.2. 
For baryonic fraction we use the values Q b = 0.05 and 0.1, 
which is in the range predicted by big bang nucleosynthe- 
sis 0.007/iT 2 <, Q b & 0.024ft- 2 (Walker et al. 1991; Copi, 
Schramm & Turner 1995). The values of the model param- 
eters are summarized in Table 1. 

2.2 Initial Power Spectrum 

Assuming the density distribution is Gaussian random, the 
initial condition is determined by the power spectrum only. 
It is given by the following functional form 

P(k,a,Q h ) =a 2 Ak n T 2 (k,a,Q h ), (1) 



© 0000 RAS, MNRAS 000, 000-000 



Baryonic fraction in the CHDM universe 3 



where a = (1 + z)^ 1 is the expansion parameter, A is a nor- 
malization constant, and k is the wave number. T(k,a,flh) 
is the transfer function describing linear changes in the per- 
turbation. We adopt the transfer function derived by Ma 
(1996). For CDM, it is given by 

ln(l + 2.34g) 



T c {k,a,Q h ) = 



2.Uq 



1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 
1 + 0.01647a; 3 ' 259/2 + 2.803 x 10~ 5 a; 3 ' 259 



1/4 



1 + 10.90a; 3 r 



1/2 



(2) 



where F = exp(-2Q b )h, q = k/Yh, Y v = a 1/2 £l h h 2 , x = 
k/Y v , xo — x(a — 1), and k is in unit of Mpc -1 . For HDM, 
it is given by 

T h (k,a,Qh) =T c (k,a,n h ) (3) 
exp(-0.0015a;') ! " 



1 - 0.121x a /' 2 + 0.102a;' - 0.0162a;' 3 / 2 + 0.00171a;' 2 



where x' = k/Y v h. We assume that CDM and baryonic 
matter have the same power spectrum. Then, the density- 
weighted power spectrum P = [fihP^ 2 + (1 — Qh)Pc^ 2 ] 2 
describes the gravitational perturbations contributed by all 
the matter components. 

The normalization constant A is related to as = a(R — 
8ft _1 Mpc, a — 1, Q h ) by the following relation 



a (R,a,Q h ) 



r — 

Jo k 



4irk 3 P(k,a,n h )W 2 (kR) 



(4) 



where W{kR) = 3[sin(fc.R) - kRcos(kR)]/(kR) :i is the top- 
hat window function. After COBE, the rms quadrupole 
Qrms-ps inferred from the cosmic microwave background 
anisotropy is often used to fix as- In the CHDM universe, 
as is approximately given by 

_ , , h 2 0.008 (1 -™ )/2 
Qis(n) 



cr 8 



0.0136 + 0.294ft - 803 + 0.109ft 2 

^(l-n^-^^.ne^fe)- - 893 

l+0.116(n h /i)-°- 893 



(5) 



where Qis(n) = Q rms _ PS (n)/18/iK (Ma 1996). We take 
Qrms-ps = 18/iK for n = 1, Q rm s-ps = 19.2^K for 
n = 0.9, and Q rms -ps = 20.5/iK for n = 0.8, respectively 
(Gorski et al. 1994; Ma 1996). 

2.3 Numerical Method 

The simulations have been performed with the cosmologi- 
cal hydrodynamic code described in Ryu et al. (1993). It is 
based on the Total Variation Diminishing (TVD) scheme, 
which is an explicit, second-order, Eulerian finite difference 
scheme (Harten 1983). Extra care was taken with the code 
in two aspects: (1) Strong shocks (with Mach number larger 
than ~ 100) were identified and handled specially in order 
to prevent the pre-shock regions from being heated unphys- 
ically. (2) Self-gravity and cosmological expansion were in- 
cluded in the way that the total energy conservation repre- 
sented by the Layzer-Irvine equation (Peebles 1980) is pre- 
served. 



Since the core radius of typical X-ray clusters is 
Si 0.5/i _1 Mpc and the cluster-cluster separation is ~ 
50/i _1 Mpc, ~ 100 3 is the absolute minimum number of 
grids required for baryonic matter in order to get statis- 
tically meaningful X-ray quantities. The simulations have 
been done in a periodic computational domain using 256 3 
grids, 128 3 CDM particles, and 2 x 128 3 HDM particles. The 
initial condition for the HDM particles has been generated 
in pairs with random and opposite thermal velocities in or- 
der to sample neutrino phase space distribution (Klypin et 
al. 1993) . The comoving size of the computational domain is 
80/i _1 Mpc, so the grid size is 0.31ft _1 Mpc. While a smaller 
size would allow us to resolve cluster structure better, a 
larger size would allow larger waves included and give us a 
larger sample of high-luminosity and high-temperature clus- 
ters. With the domain size of 80/i _1 Mpc, which has been de- 
termined by a compromise between the two considerations, 
clusters are identified comfortably and yet enough clusters 
are yielded for statistical analyses. However, with the grid 
size of 0.31/i _1 Mpc, clusters are not "fully" resolved result- 
ing in under-estimate of X-ray luminosities (see §3.2 for 
further discussion). The simulations have been started at 
z — 15 (or a — 1/16) and run to z — (or a = 1). 

In the simulations, atomic processes such as heating and 
cooling have not been included, as they will have little effect 
on the hot cluster gas discussed in this paper. 



3 RESULTS 

3.1 Power Spectrum 

We have calculated the linear power spectrum of total (dark 
and baryonic) matter for the 4 models listed in Table 1 us- 
ing the formulae in §2.2. Figure 1 shows them as well as the 
power spectrum reconstructed from galaxy and cluster sur- 
veys by Peacock & Dodds (1994). The upper panel contains 
the power spectra of the models with h/n/Qh = 0.5/0.9/0.1 
(Model CHDMc and CHDMd) and the lower panel with 
h/n/n h = 0.5/0.9/0.2 (Model CHDMa and CHDMb). They 
agree well with the Peacock & Dodds's reconstructed power 
spectrum as was already shown by Ma (1996), although 
the Peacock & Dodds's power spectrum is for the CDM 
model. Note that each panel contains a dotted curve for the 
model with = 0.05 and a solid curve for the model with 
fib = 0.1, although they are not well distinguished. In gen- 
eral, increasing Qb with fixed Oj, results in decreasing power 
at small scales, but only by a small amount. The plots indi- 
cate that the power spectrum is insensitive to the baryonic 
fraction. The power spectra of other CHDM models with 
0.5 < h < 0.8, 0.8 < n < 1, and 0.05 < fi h < 0.3 can be 
found in Choi (1996). 



3.2 Properties of X-ray Clusters 

The properties of X-ray clusters are sensitive to Qb, since 
the X-ray bremsstrahlung emission is proportional to the 
square of gas density. So they can be used to constrain Qb- 
We select two models, Model CHDMa and CHDMb with 
h/n/n h = 0.5/0.9/0.2 and Q b = 0.05 and 0.1, and study 
by numerical simulations the properties of X-ray clusters in 
those models in detail. 
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Figure 1. Linear power spectra of total (dark and baryonic) matter for Model CHDMc and CHDMd (upper panel) and CHDMa and 
CHDMb (lower panel). Each panel contains a dotted curve for the model with il;, = 0.05 (Model CHDMc and CHDMa) and a solid curve 
for the model Qj, = 0.1 (Model CHDMd and CHDMb). The reconstructed power spectrum from galaxy and cluster surveys is shown 
with filled circles (Peacock & Dodds 1994). 



X-ray clusters emit X-rays from the hot intra-cluster gas 
that fills the space between galaxies and occupies much of 
the cluster's volume. In general the core radius of typical X- 
ray clusters is about 0.5fa _1 Mpc and the total radius about 
l/i _1 Mpc, but they cannot be accurately determined from 
observations. The total (bolometric) X-ray luminosity due 
to the thermal bremsstrahlung is 



47T 



jffdr dv 



(6) 



where jff is given in units of erg cm 



-V^z^sr" 1 by 



Jff = TZ 



32e 4 



47T 3ml c 3 



7T/h>o(H) 



3kT 



1/2 



exp 



V kTJ 



x g ff (T, v) [n(HII) + n(Hell) + 4n(HeIII)] n(e). (7) 

gff(T,v) is the Gaunt factor. We assume the primordial 
abundance, 76% H and 24% He by mass, and that H and 
He are fully ionized. 

X-ray clusters in the simulations are identified as fol- 
lows. We first calculate the total X-ray luminosity due to 
the thermal bremsstrahlung in each cell. The cells with the 
total X-ray luminosity higher than 10 38 erg s _1 are selected 



as X-ray bright cells. Then, we find the local maxima by 
comparing the total X-ray luminosity of each X-ray bright 
cell with that of 124 neighboring cells, and identify them 
as the centers of the X-ray clusters. Having identified the 
X-ray cluster centers, we define the X-ray clusters in the 
whole computational domain. The total X-ray cluster vol- 
ume consists of 125 cells, the central cell and 124 cells sur- 
rounding it. For the calculations of X-ray luminosity and 
temperature, each cell is weighted equally by a weight fac- 
tor = 47r/3 (radius) 3 /(resolution) 3 /(the number of cells) 
= 4tt/3 (1.0)7(0.31) 3 /125 = 1.098, so that the total vol- 
ume of each X-ray cluster equals the volume of a sphere 
of radius l/i -1 Mpc. The weighting scheme compensates for 
the adoption of a slightly small volume by heightening the 
weight per cell. 

We note that our scheme is similar to (but not same 
as) that used by Kang et al. (1994), Cen & Ostriker (1994), 
and Bryan et al. (1994a), but different from that used by 
Bryan et al. (1994b). Our test showed that the calculated 
luminosity and temperature are insensitive to the details 
of the scheme, as already pointed by Bryan et al. (1994b). 
However, as pointed in those previous works, our luminosity 
should have been under-estimated due to insufficient reso- 
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Figure 2. Projected distribution of the X-ray clusters with L x > 10 41 erg s 1 from the region of r < lh 1 Mpc. The open circles in the 
left panels represent the X-ray clusters in CHDMa, and the filled circles in the right panels represent the X-ray clusters in CHDMb. 



lution to resolve the central density structure in clusters. 
Bryan et al. (1994b) pointed that in their simulation with 
resolution similar to ours, the under-estimation in the lu- 
minosity could be as high as a factor of 3. The importance 
of resolution was also point by Evrard, Metzler & Navarro 
(1996) in the study of X-ray clusters using the SPH code. 
On the other hand, the calculated temperature is relatively 
insensitive to resolution. 

The bright X-ray clusters with L x > 10 41 erg s _1 from 
the region of radius l/i _1 Mpc at z — are shown in Figure 
2. Left panels show the projection of the X-ray clusters gen- 
erated in CHDMa in the x — y, y — z, and z—x planes, and the 
right panels show the projection of the X-ray clusters gener- 
ated in CHDMb. The distribution pattern of the X-ray clus- 
ters is similar, but there are more X-ray clusters in CHDMb 
(with larger flj,) than in CHDMa. In CHDMa, our computa- 
tional domain with (80fe _1 Mpc) 3 volume at z — contains 
no cluster with total luminosity brighter than 10 45 erg s _1 , 
no cluster brighter than 10 44 erg s _1 , 13 clusters brighter 
than 10 43 erg s" 1 , 44 clusters brighter than 10 42 erg s _1 , and 
268 clusters brighter than 10 41 erg s _1 . In CHDMb, our com- 
putational domain at z — contains no cluster with total lu- 
minosity brighter than 10 45 erg s _1 , 2 clusters brighter than 
10 44 erg s"\ 24 clusters brighter than 10 43 erg s"\ 126 clus- 
ters brighter than 10 42 erg s _1 , and 558 clusters brighter than 
10 41 erg s _1 . 

The luminosity functions calculated with the iden- 
tified X-ray clusters are shown in Figures 3a and 3b 
for CHDMa and CHDMb. The plots show the luminos- 
ity functions in 10 40 erg s _1 < L x < 10 45 crg s _1 and 
< z < 1. There are very few high- luminosity clusters 
with L x J> 10 44 erg s~\ because our computational do- 
main is too small to contain them. Also the plots show 
the observed bolometric luminosity function {3.1^ 4 'g x 
10- 6 /i 3 Mpc- 3 ft 2 [L 4 4(&oO]" 1 }[ft 2 L44(&oO]~ 1 ' 85±0 ' 4 by Henry 
& Arnaud (1991) as the boxed area. The Henry & Arnaud's 
luminosity function agrees with the more recently considered 
ones, such as the one by Burn et al. (1996), although White, 
Efstathiou & Frenk (1993) pointed errors in the Henry & 
Arnaud's luminosity function. The observed bolometric lu- 
minosity function is more consistent with the calculated lu- 
minosity function of CHDMa than that of CHDMb. A pos- 
sible correction for the under-estimation of the X-ray lumi- 
nosity would make the discrepancy between the calculated 
luminosity function of CHDMb and the observed luminosity 
function worse. 

Also in CHDMb we find weak positive evolution until 
z ~ 0.2, and then mild negative evolution thereafter. This 
can be seen in Figure 4 which shows the number density 
evolution of X-ray clusters with L x > 10 erg s~ between 
z — and z = 1.5. Open circles are for CHDMa and filled 
circles are for CHDMb. We note that the evolution of the 
cluster abundance does not agree with the analytic predic- 
tions such as the one based on the Press-Schechter approx- 
imation (Bahcall, Fan & Cen 1997; Fan, Bahcall, & Cen 
1997), according to which the evolution is much steeper in 
high Q and low as models. 



10 3 



o 
2 



A 
J* 




10- 1 



Figure 4. Evolution of the number density of the X-ray clusters 
with L x > 10 43 erg s — 1 from the region of r < l/i _1 Mpc. The 
open circles and dotted line arc for CHDMa, and the filled circles 
and solid line are for CHDMb. 



We fit the luminosity function to an approximate 
Schechter function 



n(L)dL = n (^r) 



dL 

T* 



(8) 



using the Levenberg-Marquardt method, and determine the 
Schechter parameters (a, L* , no) as a function of redshift. 
The parameters are calculated at seven different redshifts 
and listed in Table 2. Here, L* is in unit of 10 44 erg s _1 and 
no is in unit of 10~ 6 /i 3 Mpc~ 3 . The Schechter functions with 
these parameters are plotted with dotted lines for CHDMa 
and solid lines for CHDMb in Figure 3. The value of the 
Schechter a-parameter in CHDMa, ~ 1.7, is slightly larger 
than that in CHDMb, but still smaller than the best-fit ob- 
servational values, 1.9 — 2.0 quoted by Henry (1992). How- 
ever, the Schechter a-parameters in our simulations are pri- 
marily determined by lower luminosity clusters than those 
used for the fits of observational data (Henry 1992). The 
characteristic luminosity, L* , and the total number density, 
no, of the Schechter luminosity function are determined less 
reliably, because the number of samples is too small. 

The temperature distribution function, in addition to 
the luminosity function, of X-ray clusters provides a useful 
test for the structure formation theory. It has been calcu- 
lated at four different redshifts and shown in Figures 5a and 
5b for CHDMa and CHDMb. The turnover at low X-ray 
temperature, T x , is possibly caused by our definition of min- 
imum cell luminosity to constitute X-ray bright cells. There 
are almost no high-temperature clusters with T x <; lOkeV. 
Again, this is due to the limited volume of our computational 
domain which is not big enough to contain high-temperature 
clusters. The observed temperature distribution function 
is (lStH x 10" 3 /i 3 Mpc" 3 keV" 1 )(fcT)- 4 7±0 - 5 according to 
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Figure 3. Luminosity function of the X-ray clusters at four different epochs with X-ray luminosity from the region of r < l/i _1 Mpc 
integrated over the whole frequency range. The open circles in (a) represent the luminosity function in CHDMa, and the filled circles 
in (b) represent the luminosity function in CHDMb. The dotted and solid lines are the Schechter fits, and the boxed area shows the 
observational data (Henry & Arnaud 1991). 
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Figure 5. Temperature distribution function of the X-ray clusters at four different epochs. The open circles in (a) represent the 
distribution function in CHDMa, and the filled circles in (b) represent the distribution function in CHDMb. The boxed area shows the 
observation data (Henry & Arnaud 1991). 
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Table 2. Schcchter fits for the X-ray clusters luminosity function. 
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Figure 6. Relation of T x versus L x at z = 0. The open circles 
in the upper panel are from CHDMa, and the filled circles in the 
lower panel is from CHDMb. The dotted and solid lines show 
the best fits, and the boxed area indicates the observational data 
(Henry & Arnaud 1991). 



Henry & Arnaud (1991), and shown as the boxed area. In 
the reliable range of T x , the distribution functions for the 
both models are almost identical. This is expected, since T x 
is mainly determined by the gravitational potential of total 
matter. Both the distribution functions agree well with the 
observed data. 

The relation between the cluster's temperature, T x , and 
luminosity, L x , should be less sensitive to the normalization 
of power spectrum. In Figure 6, we show the plots of T x 
versus L x at z = 0. The upper panel is for CHDMa and the 
lower panel is for CHDMb. The boxed area indicates the 
observed relation, log 10 T x (heV) = log 10 (4.2t + (0.265 ± 
0.035) log 10 (fr 2 L 4 4) according to Henry & Arnaud (1991). 
The best straight-line fits, log 10 kT x = A + B\og 10 L x , 
are shown as a dotted line for CHDMa and solid line for 
CHDMb. A = 0.672 and B = 0.413 for CHDMa, and 
A = 0.495 and B = 0.466 for CHDMb, respectively. For 
a given luminosity, the temperature is somewhat lower for 
CHDMb (A = 0.495) than for CHDMa (A = 0.672), and 
the slopes (B = 0.413 and 0.466) are somewhat steeper than 
that indicated by observations (B = 0.265±0.035). However, 
the observed slope is determined mainly by high-luminosity 
clusters. We find that in the region where comparison can 
be made, the agreement with the observations is better for 
CHDMa than CHDMb. 

For each simulated X-ray cluster, the gas mass and the 
total mass within a sphere of radius l/i _1 Mpc are calcu- 
lated at z = 0. The histogram of the ratio of the two 
masses is shown in Figure 7. The thick dotted histogram 
shows the ratio distribution for CHDMa, and the thick 
solid shows the ratio distribution for CHDMb. These his- 
tograms are arbitrarily normalized to have a peak height 
similar to the observation data, which is drawn with the 
thin solid histogram, adopted from Jones & Forman (1992). 
The more recent observed ratio of the gas to total mass 
is (Mgas/Mtot)obs = 0.1 — 0.22 for a refined sample of 
13 clusters (White & Fabian 1995), so the observed ra- 




M gas /M lol 

Figure 7. Histogram of the ratio of gas to total mass in the X- 
ray clusters with L x > 10 41 crg s _1 at z 



0. The thick dotted 
line represents the histogram from CHDMa, and the thick solid 
shows the histogram from CHDMb. The thin solid line shows the 
observational data (Jones & Forman 1992). 



tio has gone up a little. The median of the computed ra- 
tio is {M gas /M t ot)com = 0.0731 ± 0.0038 for CHDMa and 
(M gas /Mtot)co m = 0.120 ± 0.0069 for CHDMb. Improving 
observations will probably narrow the observed histogram, 
while increasing the dynamical range of numerical simula- 
tions and incorporating more realistic physics will widen the 
computed histograms. Both expected improvements should 
make the agreement better. 

Figure 8 shows the evolution of the average X-ray tem- 
peratures for clusters with L x > 10 43 erg s _1 from z — to 
z = 1.5. Open circles are for CHDMa and filled circles are 
for CHDMb. Also shown are the best-fits to the evolution 
analytically predicted by Kaiser (1986), T x oc (1 + z)^ 1 for 
n = 1. The evolutions for the both model are similar, as 
expected. But there is an overall difference that the aver- 
age temperatures are smaller in CHDMb by a factor of 0.65 
than those in CHDMa. The computed evolution of the av- 
erage temperatures agrees only marginally with analytical 
prediction by Kaiser (1986). 



4 CONCLUSION 

We have computed the linear power spectrum of total mat- 
ter for spatially flat, COBi^normalized CHDM models. We 
have found that the models with h/n/Qh = 0.5/0.9/0.1 and 
0.5/0.9/0.2 give a power spectrum in good agreement with 
the observed power spectrum. Then, through numerical sim- 
ulations that include the hydrodynamics of baryonic matter, 
as well as the particle dynamics of dark matter, we have 
computed the properties of X-ray clusters for the models 
with h/n/n h = 0.5/0.9/0.2 and Q b = 0.05 and 0.1. We have 
found that the models with Q b = 0.05 produce X-ray clus- 
ters with properties that agree well with the observed data. 
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Figure 8. Average X-ray temperature for the clusters with L x > 
10 43 erg s — 1 as a function of redshift. The open circles in the 
upper panel is from CHDMa, and the filled circles in the lower is 
from CHDMb. The dotted and solid curves show the best-fits to 
the evolution predicted by Kaiser (1986), T x oc (1 + z) ~ 1 . 



Our results have shown that the physical properties of 
X-ray clusters considered in this paper are quite sensitive 
to rij, , since the bremsstrahlung emission is proportional to 
pf. Although the number of X-ray bright clusters in the 
model with Qb = 0.1 exceeds slightly the observed number, 
we still consider the model to be marginally consistent with 
observations. However, the COBi^normalized CHDM model 
with Qb > 0.1 may be ruled out by the present work. Finally, 
we note that that our simulations are not fully converged, 
so X-ray luminosity may have been under-estimated by a 
factor of 3 or more. Hence, future simulations with higher 
resolution will put a more stringent constraint on fib- 
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